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I. INTRODUCTION 

Topological phases of condensed matter have attracted 
increasing attention in recent years. Examples include 
spin liquids, fractional quantum Hall (FQH) states, and 
topological insulators, which exhibit fractionalized exci- 
tations and statisticaii^i^i^i^i^iiiSiiaiiiii^ii^iii. A partic- 
ularly novel application has been the proposal of uti- 
lizing non-Abelian elementary excitations of topologi- 
cal phases to achieve fault-tolerant topological quan- 
tum computationiii^. To this end, Kitaev, in a semi- 
nal paper—, introduced a model of interacting spins on 
a honeycomb lattice which reduces to the problem of 
Majorana fermions coupled to a static Z2 gauge field. 
The ground state is topologically nontrivial and breaks 
time reversal (TR) symmetry, and the elementary exci- 
tations are anyonic, with non-Abelian statistics. Much 
progress has been made toward understanding the na- 
ture of this phas6a^^^^^i20^^2-23,24^ -^^ 

eralizing the model to other lattices^ '^^i^^ and to three 
dimensions^i. Very recently, a generalization to a dec- 
orated honeycomb lattice exhibiting a chiral spin liquid 
ground state was made by Yao et a^. On the experi- 
mental side, the y = \ FQH state is expected to real- 
ize non-Abelian anyons^SiSSi^. Time-reversal invariant 
topological states have also been realized in HgTe semi- 
conductor quantum wells^ and in Bii.^^Sbj^ alloys^. 

The solvability of the Kitaev model depends crucially 
on the property of the Pauli-matrices, i.e., {ai,aj} — 
2dij and (JxUyaz = i, which are the simplest example of 
a Clifford algebra. This gives rise to the constraint that 
all the above models are defined in lattices with the co- 
ordination number three, and thus most of them are on 
lattices of dimension two. Extending the Kitaev model 
to more general lattices, three dimensions, and large spin 
systems enriches this class of solvable topological models. 
These extensions naturally involve higher ranked Clif- 
ford algebras, with 2" x 2"(n > 2)-dimensional matrices, 
which can be interpreted as high spin multipole opera- 
tors. Some early work on exactly solvable models in the 
F-matrix representation of the Clifford algebra has been 
done in refs. [33,34,35]. 

In this article, we generalize the Kitaev model from 
the Pauli matrices to the Clifford algebra of F-matrices. 



For the 4x4 representation, we construct a model in a 
decorated square lattice with the coordination number 
5, which can be interpreted as a spin-| magnetic model 
with anisotropic interactions involving spin-quadrapole 
operators only. It is interesting that although each spin- 
quadrapole operator is TR invariant, the ground state 
spontaneously breaks TR symmetry. Such a state is a 
topologically nontrivial chiral spin liquid state with ex- 
tremely short ranged spin correlation functions. The 
topological excitations are expected to be non-Abelian. 
The F-matrix formalism is also convenient to define a 3D 
counterpart of the Kitaev model on the diamond lattice. 
By breaking the TR symmetry explicitly, a gapless spin 
liquid with a 3D Dirac-cone like spectrum is found. Topo- 
logical insulating states with TR symmetry also may be 
elicited on the diamond lattice. A generalization to even 
higher rank F-matrices is also discussed, wherein a topo- 
logically non-trivial spin-liquid state appears along which 
manifests a suitable defined "time-reversal-like" symme- 
try. 

II. 2D CHIRAL SPIN LIQUID WITH 
F-MATRICES 

A. Remarks on the nature of Kitaev's model 

Before elucidating the details of our model, it will prove 
useful to reflect on why the Kitaev model is equivalent to 
noninteracting fermions in a static Z2 gauge field. The 
F-matrices obeying the Clifford algebra {F",F''} = 2(5°'' 
may be represented in terms of 2n Majorana fermions 77 
and {1^"}, where a — 1, ... , 2n— 1. Then define the 2n— 1 
F-matrices F° — i-qS,"^ . The product 

A = Pip^ . . . f2"-i = ir;^^^^ • • • ^2"-! (1) 

then commutes with each of the F° and furthermore satis- 
fies A^ = (—1)""-'^, hence we can choose A = A = 
which is to say that all states in each local Hilbert space 
satisfy Aj^*) ~ Aj^*). Now consider a lattice C of coordi- 
nation number z = 2n—l in which the link lattice is itself 
z-partite. That is to say that each link can be assigned 
one of z colors, and no two links of the same color termi- 
nate in a common site. On each link (ij) of the lattice. 



then, we can write interaction terms Tf = —u-- -in-n-, 
where i and j are the termini of the hnk, and u^^ = i^f , 
where a is the color label of the link. Note that m?- = 1, 

and furthermore [WyjUfc;] = for all {ij) and {kl), and 
in particular even if these links share a common termi- 
nus. Thus, the set {w^} thus defines a configuration for 

a classical TL-i gauge field. Note that m-^- — —u^^. The 
model, is then 



D' C" 



(2) 
(3) 



where the 'color' index a is associated with each partic- 
ular link. The honeycomb lattice satisfies the above list 
of desiderata, with z = 3. Accordingly, the Kitaev model 
has interactions crfa-j on 0° links, crfuj on 120° links, 
and af o-f on 240° links. 

For n = 2, the F matrices are of course the Pauli ma- 
trices, and the set {l,(T^,(Ty ,a^} forms a basis for rank-2 
Hermitian matrices. For n = 3, in addition to the five F° 
matrices, we can define (^"2^^) = 10 additional matrices, 



(4) 



resulting in a total of 1 + 5 + 10 = 16, which of course 
forms a basis for rank-4 Hermitian matrices. For n = 4, 
we add the (^"3"^) = ^6 matrices F'^''^ = 77^"^'^ to the 
F''^ (21 total), F" (7 total), and 1, resulting in the 64 
element basis of rank-8 Hermitian matrices, etc. 

For each closed loop C on the lattice, one can define a 
Z2 flux, which is a product. 



Fc = 



(5) 



where the product is taken counterclockwise over all links 
in C If C is a self-avoiding loop of N sites (and hence N 
links), then 



- -F" 



•F 



(6) 



These fluxes are all conserved in that they commute 
with each other and with the Hamiltonian. Under a lo- 
cal gauge transformation, the Majoranas transform as 



-77j, which is equivalent to taking 



for 



each link emanating from site i. The gauge-invariant con- 
tent of the theory thus consists of the couplings {Jij} and 

the fluxes {F^,} associated with the elementary plaque- 
ttes V. For a given set of fluxes, there are many (gauge- 
equivalent) choices for the u^y By choosing a particular 
such gauge configuration, the Hamiltonian of eqn. [3] may 
be diagonalized in a particular sector specified by the Z2 
fluxes. 

It is worth emphasizing the following features of the 
as defined in eqn. [5] First, a retraced link contributes a 
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FIG. 1: (color online) The decorated square lattice (right) 
with linked diagonal bonds for the Hamiltonian of eqn. ??. 
Each unit cell (left) consists four sites (A, B, C, D), ten links, 
and six plaquettes including two squares and four triangles. 
The bonds are classified into five types as marked with labels 
a which run from 1 to 5 and 1' to 5'. 



factor of —1 to the flux, because u,, ■ u,, = — 1. This has 
consequences for combining paths. If two loops C and C 
share k links in common, thenF^c- = (-1)'=^^^^,, where 
CC is the concatenation of C and C , with the shared links 
removed. Consider, for example, the triangles ABD' and 
BC" D' in the left panel of fig. [TJ The combination of 
these triangles yields the square ABC" D' . One then has 



ABD 



BC"'D' 



'■AB"'BD'"-D'A 



BC'ii'^C'" D^DB 
(-1) 

(7) 



AB'^BC 



,,,U 



C"'D''^D'A 



= -F 



ABC'D' 



since the single link BD is held in common, but is tra- 
versed in opposite directions. 

A second point regarding the fluxes is that if n(C) 
is the number of links contained in C, then traversing C 
clockwise rather than counterclockwise results in a flux 
of F^_i — (— 1)"(''^F|^. Thus, the Z2 flux reverses sign for 
odd length loops if the loop is traversed in the opposite 
direction. For loops of even length, the flux is invariant 
under the direction of traversal. 



B. Model definition 

Our generalization of Kitaev's model involves an n = 3 
system on the decorated square lattice depicted in fig. [TJ 
The link lattice is 5-partite: each lattice node lies at the 
confluence of 5 differently colored links. We will take the 
side length of each square to be a. The lattice constant 
for the underlying Bravais lattice, i.e. the distance be- 
tween two closest A sublattice sites, is then a' — 2a. The 
model is that of eqn. [21 



n = -,hJ2 r] - J2 ^ f^ (s) 

1 links 2 links 

-^Er?r3-j,;^FfF|-j5;^F^F^ 



3 links 



4 links 



5 links 



While the structural unit cell shown in the left panel 
of fig. [H contains distinct bonds, labeled with a among 
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{1, 5} or {!',..., 5'}, and is solvable with 10 distinct 
couplings J^, we shall assume = Jy etc.. 

The 4x4 F-matrices may be explicitly taken as 



-I 

1 



, r' 



2,3.4 



a 

-CT 











.(9) 

where I and a are the 2x2 unit and Pauli matrices, 
respectively. These five F-matrices are in fact spin- 
quadruple operators for the spin-| system^. The other 
ten F-matrices, defined above in eqn. |4j contains both 
spin and spin octupole operators22i2S. For later con- 
vience, the ten F"'' matrices are also written explicitly 
here as 



t12,13,14 



-^34,42,23 



a 

a 

a 

a 



F 



25,35,45 



p51 ^ 



-a 
3 



I 

-I 



(10) 



From the six Major ana fermions, we may fashion three 
Dirac fermions, viz. 



^04 



(11) 

(12) 

C23 = w-^e) ■ (13) 

The three diagonal F-matrices can be represented as 

1 (14) 



T^4 _ n t 

J- — ^t-04''04 



Fi5 = 2c{,c 



5^15 



1 



p23 _ n t _ 1 

^ — ^''23''23 ^ 



(15) 
(16) 



The basis vectors of the physical space in terms of the 
S'^-eigenstates are: 



2) ^ 43''15''04 



n) 



-15 



1^) 



\) = 43 1^^) 



(17) 
(18) 
(19) 
(20) 



where jfi) is the reference vacuum state. In other words, 
the physical states have odd fermion occupation number 

The i?-matrix defined in Sect. 



04 '-04 "rC]^5C]^5-|-C23C23 



III CI can be represented as i? = —^^15 + ^i5)i^: 



23 



For each of the six plaquettes per unit cell, one defines 
a Z2 flux as in eqn. \5\ For example, for the ABD' triangle 
shown in fig. [1] enclosed by the a = 2, a = 3, and a — 5 
bonds, we have 



^ABD' 



{&B){&D'){^eD'eA) 



r23 7^35 "p52 
A ^ B ^ D' 



(21) 



These fluxes are all conserved in that they commute with 
each other and with the Hamiltonian. As pointed in refs. 



the flux configuration for the ground state on the 
triangular lattice is odd under TR symmetry, thus the 
ground state breaks TR symmetry and is at least doubly 
degenerate. 



C. Time reversal 

The time reversal (TR) transformation is defined as the 
product T — RC , where C is the complex conjugation op- 
erator and R is the charge conjugation matrix, satisfying 
R^ = -I R'^ = R-^ = R^ = -R. Explicitly, we can 
take R = F^F^ = Note that RT" R = -(F°)* and 

RT"-'' R= (F"'')*. Acting on the Majoranas, the complex 
conjugation operator C is defined so that 









1 ,L \ 






H 










[ ■« j 



(22) 



With these definitions, the F° operators are even under 
TR while the F"'' operators are odd. Note that RrjR — rj. 
Thus, the effect of the time reversal operation on the 
fluxes is the same as that of reversing the direction 
in which the loop is traversed, i.e. RF^^R — (— 1)"*^''-'F(^, 
where n{C) is the number of links in C. 



D. Projection onto the physical sector 

As we have seen, in terms of the Majorana fermions, 
the F-matrices are represented as 



F'^ = ^77^° 



(23) 



We further demand that any physical state Y^) must 
satisfy AJ*) = -l^-), where A, = F^FfFfFfFf, on 
each lattice site i (see eqn. [Ij. That is, each state in 
the eigenspectrum is also an eigenstate of the projector 
P — rii 5(1 ~ ^i); SO the local Hilbert space at each 
site is four-dimensional, rather than eight-dimensional. 
Since [i?, P] = 0, the two operators can be simultane- 
ously diagonalized. For any eigenstate \^) oi H in the 
extended Hilbert space (with local dimension eight), we 
have that is also an eigenstate of H , and with the 

same eigenvalue. Thus, we are free to solve the problem 
in the extended Hilbert space, paying no heed to the local 
constraints, and subsequently apply the projector P to 
each eigenstate of H if we desire the actual wavefunctions 
or correlation functions. Note also that while the projec- 
tor P does not commute with the link Z2 gauge fields 
, nevertheless [P, F^] — for any closed loop C, so the 
projector does commute with all the Z2 fluxes, which, 
aside from the couplings, constitute the gauge-invariant 
content of the Hamiltonian. 
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E. General Majorana Hamiltonian 

The general noninteracting lattice Majorana Hamil- 
tonian is written as H = ^ J2i j ^ij Vi Vj ; "^ith H = 

= -H* and {7?j,?7j} = 26^j. Let R denote a Bra- 
vais lattice site, i.e. a site on the A sublattice, and 
the index I E {!,..., 2r} labels a basis element; om- 
model has r = 2. The Majorana fermions satisfy 
{'rii{R),rii,{R')} — 2(5^^,(5;;,, and Fourier transforming 

to r]i{k) = N-^^^J2jir]i{R) e-'''^-(^+'^i\ where N is the 
number of unit cells and dj is the location in the unit cell 
of the Z'^ basis site, we arrive at the Hamiltonian 

^ = 5E^"'('^)'?/(-^)'?;'(^) ' (24) 

k 

where -ff„,(fc) = J^R^wiR - R') e'''<^'-^+'^i'-'^''> sat- 
isfies 

H^,{k) = Ht,,{k) = -H,,{-k) = -Hti,{-k) , (25) 

and rii{—k) — r]j{k). The eigenvalues of H{k) occur in 
pairs {Ej{k) , —Ej{—k)^, where j S {1, . . . , r} where 2r 
is the number of basis elements. Written in terms of 
Dirac fermions, the diagonalized Hamiltonian takes the 
form 

^ = EE(27l.7,.,-l)i?,(fc), (26) 
k j=l 

and therefore the ground state energy is 

r 

Eo = -T.T.\^M (27) 

k j=l 

and the excitation energies are t^j{k) — 2\Ej{k)\. 

F. Bulk and edge spectra 

Following the procedure outlined by Kitaev^', we repre- 
sent the Hamiltonian H in the extended Hilbert space in 
terms of free Majorana fermions hopping in the presence 
of a static Z2 gauge field, as in cqn. [31 The spectra can be 
solved in each gauge sector with a specified distribution 
of the Z2 phases Uij = ±1. 

In the extended Hilbert space, all the link phases m^- 
mutually commute with each other and with the Hamil- 
tonian. The product of the link phases around a given 
plaquette gives the Z2 flux associated with that plaque- 
tte, as in eqn. [51 It is the Z2 fluxes of all the triangular 
and square plaquettes which define the gauge-invariant 
content of the model. 

The structural unit cell, as shown in the left panel of 
fig. m consists of four sites, ten links, and six plaquettes. 
Suppose that the arrangement of fluxes has the same 




FIG. 2: (color online) When the magnetic and structural unit 
cells coincide, there five Z2 gauge degrees of freedom, which 
we take to be cti_5 as shown, with uad = etc. The other 
Uij are positive for links where the black arrow points from i 
to j. 

period as this structural cell. What of the link phases 
Mjj? It is easy to see that if the total Z2 flux of the 

structural unit cell is -1-1, i.e. if the product of the 
over the six plaquettes in the structural unit cell is +1, 
then the u^j may be chosen so as to be periodic in this 
unit cell. In such a case the magnetic unit cell coincides 
with the structural unit cell.^ If, however, the net Z2 flux 
per structural unit cell is —1, then the smallest magnetic 
unit cell (i.e. periodic arrangement of the link phases 
Mjj) necessarily comprises two structural cells. Here, we 
assume that the net flux per structural cell is -f 1, so the 
magnetic and structural unit cells coincide. Thus, there 
are five (and not six) Z2 degrees of freedom per unit cell, 
which, using the labels of fig. [2l can be taken to be the 
link phases 

Uad = ^1 ' ^Ac = ^2 . ^AB' = ^3 (28) 
^AD' = ^4 : Ugjj, = 0-5 . (29) 

The remaining five values of u^^ can then be fixed, and 
we take 

^AB — ^CB ~ ^C'B' — '^CD — ^CD' ~ ^ ■ (^0) 

The fluxes of the triangular and square plaquettes are 
then given by 

^ABCD = ^BC"'D' = ^5 

^AD'C'B' — ^3 '^i ^C'DA = ^^l'^2 

^ABD' = ^^4, ^AB'C ^ ^2^3 ■ 

Thus, there are 2^ = 32 possible distinct flux configura- 
tions which are periodic in this unit cell. Any arrange- 
ment of the Uij phases consistent with zero total flux per 



5 



cell is therefore identical or gauge-equivalent to one of 
these 32 configurations. 

We set the coupling on the horizontal links to Jh, on 
the vertical links to Jy, and on the diagonal links to Jd- 
The independent nonzero elements of H^, (fc) are then 







(31) 


HAcik) 




(32) 






(33) 




= -2i Jv cos(6'2/2) 


(34) 






(35) 


He oik) 


= 2i Jh cos(6'i/2) , 


(36) 



where we define the angles {O^.O^) — (k^a' , fc^a'). 



We have found that the ground state energy is mini- 
mized when the Z2 flux through each of the four squares 

-1, and the triangles all are 
: ±1. Under time reversal 



in the unit cell is -Fj-| 
the same value, i.e. F, 



A 



^ and 



'F^. Note that the Z2 flux 



through either square which contains a diagonal bond 
is = —F^ = — 1 if the triangles contain the same 



flux. The flux pattern with F^ 
with the choice 



F 



A 



-1 is achieved 



(c^i><^2>cr3'Cr4'f^5) = (-1,-1, +1,-1,-1) • (37) 



With this link flux assignment, the Hamiltonian matrix 
H{k) takes the form 



H{k) 



I 2i Jh cos(6ii/2) JDe-^(^i+^2)/2 -2i Jv cos(6'2/2)\ 

-2i Jh cos(6li/2) -2i Jy cos(6'2/2) Jd e*(''2-'^i)/2 

i JDe*(''i+''2)/2 2i Jv cos(6'2/2) 2i Jh cos(6ii/2) 

V 2i Jv cos(6i2/2) i JDe*(''i"''2)/2 „2i Jh cos(6ii/2) / 



(38) 



The eigenvalues of iJ(fc) are found to be 



Fi_2(fc) 
where 



-2^/41 + (J2 + j2)/(fc) ± J^^Jl + J2 5(fc) 

(39) 



+ (J2 + j2)/(fc) T J^^Jl + J2 g(fc) , 

(40) 



g{k) — cos {\Q 



72 , t2 

cos {\B^ 



(41) 
(42) 



In fig. [3] we plot the total energy per site for our model 
for all possible flux configurations for our model (32 in 
total), where we have taken Jh = Jv = Ji which we 
here and henceforth assume. We explore the properties 
of our model as a function of the dimensionless param- 
eter Jd / J. Since time-reversal has the effect of sending 
F — > —F or all odd-membered loops, every state must 
have an even-fold degeneracy. We find that the flux con- 
figurations with the two lowest lying total energies are 
each twofold degenerate, and the third-lowest lying total 
energy flux configuration is eightfold degenerate. In flg. 
|4]we show the spectra E^ik) for the lowest energy flux 
configuration. 

Our model exhibits a topological phase transition as 
increasing Jjy/ J. At Jd/J = 0, the dispersion is gapless. 



with a Dirac cone at (6*1, 6*2) = (7r,7r). As J^i/ J 7^ 0, this 
spectrum becomes massive with a gap A(7r,7r) — 2Jd. 
This is a topologically nontrivial phase, characterized by 
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FIG. 3: (color online) Total energies (per site) versus a = 
tan~^(JD/J) for our decorated square lattice model Curves 
for all possible flux conflgurations are shown. The flux con- 
figurations for the two lowest-lying total energy states are 
twofold degenerate owing to time-reversal, which reverses the 
flux in the odd-membered loops. The third-lowest lying total 
energy state has an eightfold degenerate flux conflguration. 
At a = i?! the horizontal and vertical bond strengths vanish 
and the system becomes a set of disconnected dimers. 



FIG. 4: The band structure of Eq. |38] at Jd/J > 0. (A) 
Jd/J = 1 with a massive Dirac spectrum at with 
the gap value 2Jd; (B) Jd/J = 2-\/2 where a gapless Dirac 
cone appears; (C) Jd/J — 3.5. The gap value at (0,0) is 
2|Jd — 2y/2J\. A topological phase transition occurs from a 
topological nontrivial phase at Jd/J < 2-\/2 to a trivial phase 
at Jd/J > 2^/2. 

non-vanishing Chern numbers, with concomitant gapless 
edge modes inside the gap between bands 2 and 3 in 
a sample with an open boundary, as depicted in Fig. 
[51JA). As Jd/J increases, the second and third bands are 
pushed toward each other at the Brillouin zone center. 
These bands eventually touch when J^/J = 2-\/2, form- 
ing there a new gapless Dirac cone. For Jd/J > 2-\/2, the 
system is in a topologically trivial phase in which the edge 
states no longer exhibit a spectral flow, and they remain 
confined within the band gap, as shown in Fig. [5l^B). In 
this case, the edge states are no longer topologically pro- 
tected, and they are sensitive to local perturbations, and 
indeed can be eliminated by sufficiently strong such per- 
turbations. The gap located at (6*1, 6*2) = (O7O) found 
to be A(0,0) = 2|Jd - 2V2J|. 

G. Spin correlation functions 

As is the case with Kitaev's model^^, the spin correla- 
tion function in our ground state is short-ranged due to 
the conserved flux of each plaquette. We take the ground 
state wavefunction |5'q) = P]"^) where j^*) is the un- 
projected state, with each taking a value ±1. Then 

(^'qIT^TJI^q) = (*|Pr?ir?jCf^j|*). Unless site i and j 




t I . I t I , 

-1.0 -0.5 0.0 0.5 1.0 

e/7t 



FIG. 5; A) The gapless edge modes of ?7-fermions across the 
band gap appear in open boundary samples in the topolog- 
ically nontrivial phase at Jd/J < 2\/2. B) The edge modes 
become nontopological at Jd/J > 2\/2. 

are linked, and unless a = b appropriate to the bond (ij), 
it is always possible to find a loop whose flux is flipped by 
£,i or ^j. In that case, since rjirjj and P do not change the 
loop flux, the above expectation value vanishes. A simi- 
lar reasoning shows that any two point correlation of the 
type of (*c|rf F^'^I^'q) or (*(3|rf FJ'^I^q) vanishes as 
well. 

Generally speaking, a non-local spin correlation func- 
tion is finite only if, in the Majorana representation of its 
operators, the ^'s can be expressed as a product of the 
Z2 gauge phases defined on different links, i.e. 

n^"^"=n(-o (43) 

(ij) {ij) 

Similarly to the work of refi^, it is natural to expect that 
each topological excitation of a Z2-vortex traps an un- 
paired Majorana mode. These vortex excitations will 
exhibit non-Abelian statistics. We are now performing 
numerical calculations to confirm this prediction. 
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III. AN EXTENSION TO 8 x 8 T-MATRICES 

The above procedure is readily extended to even higher 
rank Chfford algebras, such as the n = 4 case with eight 
Majorana fermions and seven anticommuting 8 x 8 F ma- 
trices. We choose the Majoranas to be {??'^,^"} where 
r = 1,2,3 and a = 1,2,3,4,5. We write T°- = irj^S,"- as 
before, and we define F^" = irf'^"' . The Hamiltonian is 



'n = -Y.Ja E (r^r,"-rfFf) (44) 

a links 

« {tiH - -nW,) ■ (45) 



Physical states are projected, at each site, onto eigen- 
states of the operator A = wy^ry^r;'^^^^^^'^^^^^ with eigen- 
value +i. 

We obtain the same behavior of 77^^^ modes as before, 
and a flat band of zero energy mode of rf" . Within a fixed 
gauge choice of ity, Ti, is invariant under a suitably de- 
fined TR reversal-like transformation of as Tr}iT~^ = 772, 
Tr]2T~'^ = —771, which is not related to the physical TR 
transformation. In the topologically non-trivial phase, 
the 7;^'^ edge modes have opposite chirality, which is 
robust in the absence of a perturbation of the type 
F6 = iri^rf, which breaks the symmetry under T. This 
behavior is similar to that of the ^He-B phase in two 
dimensions, which has been recently identified as a TR- 
invariant topological superconductor—. 



IV. r-MATRIX MODELS IN THE 3D DIAMOND 
LATTICE 

The F-matrix analogy of the Kitaev model is also ex- 
tendable to the three-dimensional (3D) systems. We will 
consider here a diamond lattice, which is fourfold coordi- 
nated, as illustrated in Fig. [51 In the following, we will 
first describe a model with a 3D Dirac cone and broken 
TR symmetry, and then another one exhibiting 3D topo- 
logical insulating states with TR symmetry maintained. 



A. 3D Spin liquid with Dirac cone excitations 

We first consider an ti = 3 model (4x4 F-matrices) in 
the diamond lattice with explicity time reversal symme- 
try breaking. Recall that the diamond lattice is bipartite, 
consisting of two FCC sublattices. Each A sublattice site 
is located in the center of a tetrahedron of B sites, and 
vice versa, as shown in fig. [S) We define the unit vectors 



V3 

1 
73 



(1,-1,1) 



;)3(1'1'-1) (46) 



which point from a given A sublattice site to its four B 
sublattice neighbors. The basis vectors for the underlying 




FIG. 6: The 3D diamond lattice with the nearest separation 
a, and two sublattices A (filled circles) and B (hollow circles). 
The unit vectors 61,2,3,4 are defined from each A site to its 
four B neighbors (see text). 



FCC Bravais lattice, which can be taken to be the A 
sublattice itself, are then 



64 = ^(0,1,1) 
64 = ^(1,0,1) 



(47) 



^3 = 63 - 64 = 73(1,1,0) • 



A general A sublattice site lies at i? = Tij^a^ + rijO 



2^2 

TT-gCg The two-element diamond unit cell may be taken 
to consist of the A site at R and the B site at i2+ 64. It 
is also useful to define the null vector 04 = 64 — 64 = 0. 

The Brillouin zone of the fcc-Bravais lattice is a do- 
decahedron. The elementary reciprocal lattice vectors, 
which form a basis for a BCC lattice, are 



61 = ^71 (-1,1,1) 



62 = ^7r (1,-1,1) 
^3 



(48) 



^3 = ^71(1,1,-1) , 



and satisfy a- • b- = 2nS-j, where i,j e {1,2,3}. Any 
vector in fe-space may be decomposed in to components 
in this basis, viz. 



We then have k ■ = ^'a, with 9^ — 0. 

We begin with the Heisenberg-like Hamiltonian, 

4 

'><0 = -J2jaJ2^R^R+e^ (50) 
a=l R 

4 

= E E "^-^ ^i? • * '^ii ^R+K ■ 

R Q=l 
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FIG. 7: A non-coplanar hexagonal ring within the diamond 
lattice. 



We use the tilde to distinguish operators which reside 
on the B sublattice from those which reside on the A 
sublattice. Here we have defined 

^R-^eutR+e^ (52) 

as the Z2 gauge field on the link between R and R + e^. 

Within the diamond lattice, one can identify four 
classes of hexagonal loops. Starting at any A sub- 
lattice site, we can move in a cyclical six-site path 
by traversing consecutively the displacement vectors 
{e^ , — , , —e^ , , — e^}. We label this hexagon by 
the indices (abc) . Without loss of generality we can as- 
sume a < h < c, since any permutation of these indices 
results in the same loop, traversed in either the same or 
the opposite sense, depending on whether the permuta- 
tion is even or odd, respectively. According to Lieb's 
theorem^S, the gauge flux in the ground state is —1 in 
each such loop, thus we can set Ujj = 1 on each link, for 
each a € {1,2,3,4}, because a circuit around a six-site 
loop involves three AB links, on which the Z2 gauge field 
is = +1, and three BA links, on which u^^ = — 1. 
Consider now 

^R^R+e TiJ+e -e, =^"it"We -e^' ^"^R^R+e -e, 

(53) 

and 

^R+ejR f We, = ""fi < ' ^Vr+s^ VR+e, ■ (54) 

These operators both break time-reversal symmetry, ow- 
ing to the presence of the F""^ and factors. However, 
owing to the commuting nature of the Z2 link variables 
Uj^, we may add terms such as these to Hq and still 
preserve the key feature that the Hamiltonian describes 
Majorana fermions 77 and fj hopping in the presence of a 
static Z2 gauge field. 

To recover the point group symmetry of the underlying 
lattice, we sum over contributions on each loop (abc) . We 




FIG. 8: First brillouin zone for the FCC structure, 
define the operators 

-irab _ -pa pah pfc /rr\ 
^R - R'- R+e^'^ R+e^~e^ K'^^J 

yR=^R+ejRi'R+e, ■ (56) 

The time-reversal violating term in our Hamiltonian is 
then written as 

'^i-T.Y.i^atV^' + kbV^') ■ (57) 

R a<b 

We reiterate that this model is also exactly solvable due 
to the commutativity of the link fluxes u^, and we as- 
sume = 1 for all a and R. 

Transforming to k space, we have the Hamiltonian ma- 
trix 

where, after performing a unitary transformation to re- 
move a phase ie^*'^*i+^2+^3)/4 from the diagonal terms, 
we have 

4 

u;ik)^J2^abMda-db) (59) 

a<b 
4 

^ik) = Y,h,,sm{e,-e,) (60) 

a<b 
4 

A(fc) = z^J,e»^" . (61) 

a=l 

The eigenvalues are E±{k) = ±y/uj^{k) + |A(fc)|2. 
One finds that E^{k) both vanish when A; lies at one of 
the three inequivalent X points, which lie at the centers 
of the square faces of the first Brillouin zone, at locations 
^{b2 + b^), ^{b^+b^), and ^(bj-l-bj)- Expanding about 
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the last of these, we write 9^ = vr + -01 , ^'2 = + ^^2 ' ^^"^ 
^3 — ^'^d assuming J^ = J and /i^^ = /i^^, wc find to 
lowest order in 2 3 that 



(/J 



A(fc)= J(^l+1^2-'/'3) • 

If we write k — ^{b^ + ^2) + Q, then we have 



9x 

1z 



^(-Vi+V'2+^3) 



V3 

4 



(^1 - V'2 
(l/'l + ■02 



V^3) 
V^3) 



(62) 



(63) 
(64) 
(65) 



Thus, the term w(fc) = u!{k) can be written as hk ■ g^, 
where is a unit vector lying in the (x, y) plane, and h 
is some combination of the ft.^^. The spectrum is there- 
fore that of a deformed Dirac cone, linear in the two 
directions g^^ and 2;, and quadratic in the third. These 
deformed Dirac cones can be made gapped by introducing 
anisotropy in the J-terms, say, J4 ^ Ji,2,3 = J ■ The sys- 
tem may then become a 3D topological insulator with TR 
symmetry breaking by developing chiral surface states. 

At this point we must ask: whither the Majoranas? 
Indeed, for the Hamiltonian Ti. = Hq + Hi, the Ma- 
joranas form a flat band at zero energy. To provide a 
dispersion for the branch, we define a new Hamilto- 
nian which we form from W with the replacements 
pa ^ p5a^ Equivalently, we can skip to the final form 
of Ti. in terms of rj-^ and and replace them with 

and Cij+e ' respectively. Combining these individual 



Hamiltonians to form H 
H{k) of the form 



W-n^ results in a matrix 



H{k) = 




(66) 



/ uo[k) 

-w(fc) 
A*(fc) 

V -A*(fc) 

w(fc)7''-RcA(fc)7i'* 



where the row and column indices range over 77, 77, 
and consecutively. The 7- matrices here are used as 
a basis for 4x4 Hermitian matrices, which are not to 
be confused with the spin and spin-multipole operators 
denoted as F-matrices. Nevertheless we still choose 7° 
and 7"'' taking the same values as the matrix forms of 
and F"'' defined in eqns. [51 and fTUl 

Again the spectra of Eq. [SS] exhibits gapless Dirac 
cones at the three X-points. This gapless excita- 
tion is robust if there no mixing terms between 77 and 
fermions, which would appear as terms involving 
{7^; 7"^ I T^'' 1 7'^^}; corresponding to same sublattice hop- 
ping, or {7^^, 7^^, 7^^, 7"^^}, corresponding to alternate 
sublattice hopping. The introduction of such couplings 
can produce a gap in the spectrum and give rise to spin 
liquid states. 



B. Time reversal symmetry properties 

The Hamiltonian matrix H{k) is Hermitian, and can 
therefore be expanded as 



H{k) = J2K{k)r + Y.^abik)' 



,ab 



(67) 



a<b 



where the couplings and A^^^ are real functions of their 
arguments. Due to the conditions of eqn. [251 we see that 
these couplings come in two classes. The coefficients of 
purely real, symmetric 7-matrices must satisfy A(A;) — 
— A(— fc); we call this class odd, or — . The coefficients 
of the purely imaginary, antisymmetric 7-matrices must 
satisfiy A(A;) = A(— A;); we call this class even, or -|-. Thus, 
we have 

class — : I , A2 , A4 , A5 , , A]^4 , A23 , A34 , , A35 
class + : A]^ , A3 , Xi^ , A24 , A25 , A45 . (68) 

The periodicity under translations k ^ k + G through 
a reciprocal lattice vector then requires A(G/2) — 
X{—G/2) — for the odd class. For three-dimensional 
systems, there are eight wavectors within the first Bril- 
louin zone which satisfy this condition, i.e. k — i (71,161 -|- 
71262 + ^^3^3) for rti = or 1, which are identified as the 
zone center F, the four inequivalent L-points, and the 
three inequivalent X-points (see fig. [5]). 

If we consider (77 , ^^) as a Kramers doublet of pseudo- 
spin up and down, we can define a TR-like anti-unitary 
transformation T as 



(69) 



where C is the complex-conjugation operation. Under 
this operation, we have 



T^^T-^ = -ri. (70) 



The 7-matrices then divide into two classes under time 
reversal: even (TjT^^ = 7) or odd {TjT^^ — —7). 
We find 

r-even : 1 , 7^ 7!^, 7'^ 7'', 7^' (71) 

T-odd :7^7^7^7^7^^7",7^^7^^7^^7'^ 

Since T also reverses the direction of k, sending k ^ —k, 
time reversal symmetry requires that the coefficients of 
the T-even 7-matrices satisfy A(fc) = A(— fc), while the 
coefficients of the T-odd 7-matrices must satisfy A(fe) = 
— A(— fc). Taking into account the division into even and 
odd classes, we find that time reversal symmetry, i.e. 
TTiT~^ — H, requires the vanishing of the following co- 
efficients: 



Ai — A3 



A 



13 



Al5 



A, 



A 



35 







(72) 



(forbidden by T-symmetry) 
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7 


class 


T 


V 


VT 


7 


class 


T 


V 


VT 


I 


- 


+ 


+ 


+ 


7" 


- 


- 


- 


+ 


7^ 


+ 


- 


+ 


- 




- 


+ 


- 


- 


7^ 


- 


- 


+ 


- 




- 


- 


+ 


- 


7^ 


+ 


- 


+ 


- 


7^^ 


+ 


- 


- 


+ 


7' 








+ 


7'' 


+ 


+ 






7^ 




+ 






7^^ 








+ 


7^^ 






+ 




-y35 




+ 






7^^ 


+ 




+ 




45 

7 


+ 


+ 


+ 


+ 



TABLE I; Symmetry properties of the 7-matrices. 



We may also define a parity operator V , as 



7^ = 7*^1 . 



(73) 



where T is the lattice inversion operator which inverts 
the coordinates relative to the position r — ^64. Under 
V we then have the classification 



P-even 
-P-odd 



I,7^7^7^7^^7^^7^^7'' 
7^7^7^^7^^7^^7^^7^^7'^ 



Parity also acts on crystal momentum, sending k 
As a result, 

^2 — ^2 — \i — ^2b — ^^34 = 
(forbidden by T'-symmetry) . 



(74) 



-k. 



(75) 



Note that the product VT does not reverse fc, and is 
given by 



VT = Y'" CI 



(76) 



A summary of the symmetry properties of the different 
7-matrices is provided in table [H 



C. 3D topological states with TR symmetry 

Recently, topological insulating states in 
three dimensions have attracted a great deal of 
attentio n^^i'^-^i'^^i^'^i'^'^i'^^i'^'' . Below we will show by adding 
the hybridization between 77 and ^ fermions, we can 
arrive at the gapped spin liquid states which can be 
topologically nontrivial. 

We begin with the additional hybridization term, 



R I 

4 

+ 2^ Sab i i fl+e/ fl-f 



(77) 



9ab ^ R+e^ ^ il ^ fl-f 



Note that 



^ R^ R+e^^ R+e-e, 



^ R+e„ ^ ^ R+e, — ^R^R *^fl+e„ 



^R'^R+e-e, ' ^Vr^R 



(78) 
(79) 



Assuming once again that Uj^ — 1 for all R and a, and 
further taking g^f^ = g^^^, this leads to the matrix Hamil- 
tonian 



H{k) = 



( uj{k) I3{k) A(fc) \ 
(3*ik) -w(fc) -A(fc) 
A*(fc) -w(fc) I3{k) 
V -A*(fe) f3*{k) oj{k) J 
uj{k) 7* + Re A(fe) 7I'' - Im A(fc) 7^^ 

+ Re/3(fc) 73^ + Im/3(fe) 7^4 



(80) 



where 



(81) 



a.b 



The spectra for each momentum k is doubly degenerate, 
and the two degenerate energy levels are found to be 



E^{k) = ±Vt^2(fc) + |A(fe)|2 + |/3(fe)|2 



(82) 



In the presence of the T symmetry, Im/3(fe) — 0. This 
can be achieved if we take m = and the coupling matrix 
g to be antisymmetric, i.e. g^^^ = — .9f,^- The system still 
possesses the additional parity symmetry of 7^ = 7^^ I. 
As a specific example, consider the case 



^13 ~ ^2'i 



h 




and 



5l2 = .913 = 523 = 

5i4 ^ .924 = 534 = 5 : 



(83) 



(84) 



where the elements of 5^^^ below the diagonal follow from 
antisymmetry. Then 



uj{k) ^ h sm{0-^ -&2)+ sin(6'i - ^'3) + sin(6'2 - 613) 



f3{k) = g sin 6^ + sin 6*3 + sin 6^ 



(85) 
(86) 



Note that uj{k), ReA(fc), and Re/3(A;) all vanish at 
the eight wavevectors fc — ^(ni 61 + 712 62 + "■3 ^3) for 
rii = or 1, which includes the zone center F, the four 
L-points, and the three X-points (which are the three 
Dirac points). If Ji,2,3,4 = J, the system remains gap- 
less at three AT-points because |Im A(fe)| = 0. By tuning 
J4 7^ ^^1,2,3 — J, the system can be made gapped, where 
the gap at the X-points is |Im A(fe)| = | J4 — J\. This 
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situation is the same as the Hamiltonian of the 3D topo- 
logical insulator studied by Fu and Kano^. Following 
their reasoning, whether or not the insulating states are 
topological can be inferred from the parity eigenvalues of 
7'*^ of the occupied states. The system is topologically 
non-trivial for J4 > J, in which case it exhibit surface 
modes with an odd number of Dirac cones in open sur- 
faces. 

On the other hand, if we relax the requirement of both 
T and V symmetries, and only keep the combined sym- 
metry of VT . Then all five coefficients are allowed in 
eqn. [501 and each energy level remains doubly degener- 
ate. However, the analysis in ref. |4l| no longer applies. 
We expect then a more diverse set of topological insu- 
lators, the study of which will be deferred to a future 
investigation. 

V. CONCLUSION 

In summary, we have generalized the Kitaev model 
from the Pauli matrices to the Clifford algebra of F- 
matrices. This enriches the physics of topological states. 



including the 2D chiral spin liquids with non-trivial topo- 
logical structure, as well as that of topological spin liquids 
with time reversal like symmetries. Possible topological 
insulating states on the 3D diamond lattice were also dis- 
cussed. 
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